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Mixed valence spinels provide a fertile playground for the interplay between charge, spin, and 
orbital degrees of freedom in strongly correlated electrons on a geometrically frustrated lattice. 
Among them, AIV2O4 and LiV204 exhibit contrasting and puzzling behavior: self-organization of 
seven-site clusters and heavy fermion behavior. We theoretically perform a comparative study of 
charge-spin-orbital fluctuations in these two compounds, on the basis of the multiband Hubbard 
models constructed by using the maximally-localized Wannier functions obtained from the ab initio 
band calculations. Performing the eigenmode analysis of the generalized susceptibility, we find that, 
in AIV2O4, the relevant fluctuation appears in the charge sector in cr-bonding type orbitals. In 
contrast, in LiV 204 , optical-type spin fluctuations in the aig orbital are enhanced at an incommen¬ 
surate wave number at low temperature. Implications from the comparative study are discussed for 
the contrasting behavior, including the metal-insulator transition under pressure in LiV 2 04 . 

PACS numbers: 71.27.-|-a, 71.20.Be, 71.15.Mb, 75.25.Dk 


I. INTRODUCTION 

Multiple degrees of freedom of electrons in solids, i.e., 
charge, spin, and orbital, have been one of the central is¬ 
sues in strongly correlated electron systems. Their inter¬ 
play is a source of various fascinating properties in mag¬ 
netism, transport, and optics, which, in turn, provides 
the possibility of controlling the diverse quantum many- 
body phenomenapi. In particular, in the systems whose 
lattice structures are geometrically frustrated, such inter¬ 
play becomes more conspicuous; charge-spin-orbital en¬ 
tangled fluctuations are promoted by keen competition 
between many different quantum states, and result in 
a broader range of peculiar behavior uniquely found in 
frustrated strongly correlated electrons^. 

Spinels, one of the primary minerals, provide a rich 
playground for studying such unique properties^. Among 
them, a series of vanadium spinel oxides with mixed va¬ 
lence are of particular interest, as they exhibit a wide 
range of peculiar properties, from metal-insulator tran¬ 
sition to superconductivity. In the present study, we are 
particularly interested in two cousins in the vanadium 
spinel oxides, AIV2O4 and LiV204. In both compounds, 
V cations comprise a pyrochlore lattice structure with 
strong geometrical frustration, and the threefold t2g or¬ 
bitals are partially occupied by a half-integer number of 
id electrons: (3d)^ ® in AIV2O4 and (3d)^-^ in LiV204. 
Despite the common aspects, two compounds show con¬ 
trasting behavior, as described below. 

AIV2O4 exhibits a structural phase transition at Tc = 
700 K^. Below Tc, the pyrochlore lattice of V cations is 
distorted so as to self-organize seven-site clusters called 
“heptamers”— . The low-temperature (T) insulating state 
is interpreted by coexistence of a spin-singlet state in 
each heptamer and nearly free magnetic moments at the 
rest V sites^i^. However, it remains unclear how the 


self-organization takes place from the interplay between 
charge, spin, and orbital degrees of freedom and the cou¬ 
pling to lattice distortions. 

In contrast, LiV204 remains metallic down to the low¬ 
est T without showing any symmetry breaking^. The 
compound, however, exhibits peculiar heavy fermion be¬ 
havior at low T, which is unexpected in 3d electron sys- 
tems^i^. Despite a number of intensive studiesi2r21, the 
origin of this heavy fermion behavior remains elusive. 
Interestingly, LiV204 exhibits a metal-insulator transi¬ 
tion in an external pressurei^^— . Although a structural 
change similar to AIV2O4 was suggested at this pressure- 
induced transitional, the relation between the two com¬ 
pounds and the nature of the transition have also re¬ 
mained as an unsolved issue. 

In this study, for understanding of the contrasting be¬ 
havior in these cousin compounds, we investigate the in¬ 
stability from the high-T paramagnetic phase by analyz¬ 
ing fluctuations in the charge, spin, and orbital degrees 
of freedom. For a multiband Hubbard model constructed 
on the basis of the ah initio band calculations and the 
analysis by means of the maximally-localized Wannier 
functions (MLWF), we compute the generalized suscep¬ 
tibility, which describes the fluctuations in the charge, 
spin, and orbital channels, at the level of the random 
phase approximation (RPA). We clarify the nature of 
relevant fluctuations by analyzing the eigenmodes of the 
generalized susceptibility. We find that Coulomb interac¬ 
tions enhance dominantly tr-type charge fluctuations in 
AIV2O4, while optical-type spin fluctuations in LiV2 04. 
We also discuss the pressure effect for adding theoreti¬ 
cal inputs to the mechanism of pressure-induced metal- 
insulator transition in LiV204. 

The organization of this paper is as follows. In Sec. m 
we introduce the model and method. After introducing 
the multiband Hubbard models constructed from ML- 
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WFs obtained by ab initio band calculations in Sec. Ill Al 
we present the definition of the generalized susceptibility 
calculated by RPA in Sec. IIIBI We describe how to ana¬ 
lyze the eigenmodes of the susceptibility and classify the 
fluctuations in Sec. Ill Cl In Sec. IIIII we show the results 
on the electronic structure and fluctuations enhanced by 
electron correlations for AIV2O4 and LiV204. After pre¬ 
senting the band structures and the tight-binding param¬ 
eters in Sec. IIII Al we show the results by the eigenmode 
analyses of the generalized susceptibilities in Sec. IIIIBI 
The results are discussed in detail in Sec. El Section El 
is devoted to summary. 

II. MODEL AND METHOD 
A. Multiband Hubbard model 

In order to estimate realistic model parameters for each 
compound, we start with ab initio calculations based 
on the density functional theory with the local den¬ 
sity approximation (LDA) ^^'^^ . In the LDA calculation, 
we use a fully-relativistic first-principles computational 
code, QMAS^i. See Ref. for the details of ab initio 
calculations. 

We construct the multiband Hubbard models by us¬ 
ing the MLWFs, which are obtained from the LDA band 
structures. In the present study for AIV2O4 and LiV204, 
t2g bands are energetically separated from others (see 
Fig. HD, and hence, we adopt the models for electrons 
in the t2g orbitals, similarly to our previous study for 
Cd20s207i2&. The Hamiltonian consists of two parts, the 
one-body part Hq and two-body interaction part Hi as 

H = Ho+Hi. ( 1 ) 

The one-body part consists of three terms: 

Ho=Ht+Ht,ig + Hsoi, ( 2 ) 

where Ht denotes the kinetic energy of electrons, iLtrig 
the trigonal crystal field splitting, and iLsoi the rela¬ 
tivistic spin-orbit interaction. Meanwhile, the two-body 
interaction part consists of two terms: the Coulomb in¬ 
teractions acting for electrons at the same atomic site 
and those between different sites, as 

_ ^on—site _j_ ^inter—site 

We estimate the one-body parameters in Hq by using 
the MLWFs as follows. The kinetic term Ht in Eq. @ is 
given in the form 

Ht= Y. YkX'[{R + r,)-{B! + rp,)\ 
R,p.R',p' C.C' 

where i?, C,, cr, and p denote the unit cell, orbital, spin, 
and sublattice indices, respectively; Vp denotes the posi¬ 
tion vector of the sublattice p in the unit cell, denotes 


the creation operator for the MLWF. The transfer in¬ 
tegral is given by the overlap integral between the 
MLWFs. On the other hand, the trigonal crystal field for 
the t2g manifold takes the form of 

/O 1 1\ 

iJtrig = Atrig X 1 0 1 , (5) 

Vl 1 Q) 

where the basis is taken as dyz, and dxy Besides, 
the relativistic spin-orbit interaction is given as 

^ IQ —i<Jz icTy \ 

Hsoi = Asoi^ ■ s = X \ icTz 0 -iax , (6) 

\-i(Ty i<Jx 0 j 

where the basis is taken as {dyz t)j {dyz 4-)! {dzx t)j 
{dzx i), {dxy t), and {dxy i); cr^, and cr^ denote the 
Pauli matrices for spin indices. To estimate the coupling 
constants, we use the following relations 

Atng = J dri;l^{r)H^^^{r)ij,;,^{r) {( ^ C), ( 7 ) 

= j dr^*xy^{r)H^^^{r)^yzi{r), (8) 

where ’tp(^^{r) is the MLWF and H^^^{r) is the Hamil¬ 

tonian in the Kohn-Sham equation. 

Next, we introduce the two-body interaction parts. 
The on-site part of the electron-electron interaction 
jjon-site gji Jg given by 

TTOn—site ^ \ ^ t-j 

til ~ 2 2 -^ '-Ja/ 3 -a'l 3 ' 

aj 3 ‘^a' 13 ' 

^ 2 ] ^aR+r„^ 8 R+r„^P'R+r„^a'R+r„y ( 9 ) 

Rp 

where a denotes the orbital and spin indices as a = 
(CajCa)- Assuming the rotational symmetry of the 
Coulomb interaction, we take the coupling constant in 
the form of 

Uqp.,c^' p.' = {U — 2Jii)S(^(^iSppi + Jii{S(^pd(^'p> + S(^pi6p(^'), 

( 10 ) 

where U and Jr denote the Coulomb interaction between 
the same orbitals and the Hund’s-rule coupling, respec¬ 
tively, and Sa/3 = denotes the Kronecker delta. 

In addition to the on-site part, we introduce the inter¬ 
site part of the electron interaction, ^™ter-site^ Although 
the Coulomb interaction is long-ranged, for simplicity, we 
take into account only the dominant component between 
nearest-neighbor sites: 

^inter-site ^ ^ ^ ^ (F ^ 

(r,r') 

where Ur = the local electron density at 

site r, and ^ denotes the vector connecting the nearest- 
neighbor sites. 
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The two-body interactions are incorporated at the level 
of RPA, as described in the next section. In the RPA 
calculations, we treat the coupling constants, C/, Jh, and 
V as parameters. 


B. Generalized susceptibility 


On the basis of the one-body Hamiltonian, Hq, we 
first calculate the static bare susceptibility XaPp-a'P'p' (q)- 
The definition is given in the form 


( 12 ) 

where G^i^p/.^pik) is the non-interacting Green function. 
k is the four dimensional wave vector k = {k,ujk) where k 
and LOk denotes the wave number vector and Matsubara 
frequency, respectively. Here, the matrix is defined 
for the orbital and spin indices, a and fi [a = (Caj^Ta); 
(a and da denote the orbital and spin, respectively], and 
the sublattice indices, p. 

Next, we calculate the generalized susceptibility by in¬ 
cluding the effect of electron correlations in a perturba¬ 
tive way at the level of RPA. For the on-site interaction 
^on-site gq RPA vertex function is defined 

by 


where the first and second terms correspond to the so- 
called bubble and ladder contributions, respectively. The 
Dyson equation is written in the form of 


= Xi-P (q) + E ( 14 ) 

/i /o 


where I denotes the set of {a,l 3 ,p). This gives in 

the matrix form as 


RPA 


X 


{q) = 


I-X^°HQ)kl] ( 15 ) 


where I denotes the unit matrix. 

On the other hand, for the inter-site interaction in 
Eq. (fTTI) , the mode coupling appears in the ladder contri¬ 
bution, which is not easy to handle in the perturbation. 
Hence, we omit the ladder contribution in when in¬ 
cluding the inter-site electron interaction. Consequently, 
y^PA for this case is given in the same form as Eq. m 
while replacing Eq. (ini) by 

^C(i/3iPi;/32Q2P2 ( 9 ) = l^C(i/3i;/32Q;2'^plP2 “f lpiP2 (9 )^Q;iQ;2^/Ji/32 I 

( 16 ) 

where lpip2(q) is the coupling constant of the inter-site 
interaction in the Fourier transformed representation. 


We note that the combined method of LDA and RPA 
potentially has so-called double-counting problem of elec¬ 
tron correlations. In the present study, however, we do 
not introduce the double-counting correction for the fol¬ 
lowing reason. In similar theoretical schemes, such as 
LDA-I-DMFT— , the double-counting correction is intro¬ 
duced for adjusting the energy gap between d and p lev¬ 
els when the electron correlations are taken into account 
only for d electrons in the effective model. Meanwhile, 
when the model deduced from LDA includes only the 
d levels, the double-counting correction is not relevant. 
In the present calculations, we construct the multiband 
Hubbard model only for the d orbitals, thus we do not 
consider the double-counting correction. 


C. Eigenmode analysis 

From the generalized susceptibility y^^^, we can ex¬ 
tract the information of fluctuations in the multiple de¬ 
grees of freedom as follows. According to the fluctuation 
dissipation theorem, the generalized susceptibility satis¬ 
fies the following relation 

^{^apq^Ppq} ^ ^ XaPp\a'P'p' p' p'{ q) ~ (f^) 

a'P'p' 

Here, the Fourier transform of the creation operator is 
given by cj,= (I/D) Y,r where D de¬ 

notes the system size; happ denotes a generalized ex¬ 
ternal field conjugate to c^apq^ppq^ 5 {A) represents 
the difference of the thermal average of A from that in 
the absence of the external field. If the external field 
is parallel to the Kth eigenvector e^^p of y^^^, namely, 
k-appi^) ^appi^)’ blq- (fT 71 ) is rewritten into the form of 

^icipqCppq) = x'^iq)h'^l^p{q), ( 18 ) 

where x'^{q) is the corresponding eigenvalue. This in¬ 
dicates that S{cl^pqCpp^) also becomes parallel to 
Thus, the eigenvector is regarded as the direction of the 
fluctuation d{cl^pgCpp^), and the eigenvalue x'^{q) corre¬ 
sponds to the amplitude of the fluctuation. Thus, the 
analysis of the eigenvalues and eigenvectors provides the 
information of fluctuations. 

Specifically, the charge fluctuation Snp{q) at the sub¬ 
lattice p is defined by 

Supiq) = ^ ^ E E ®Sap(9)' (19) 

\ (A / OC 

Snp{q) can be decomposed into the orbital components 
as 


Snp{q) =^Sn(;p{q), (20) 

C 
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where C denotes the orbital index. Here, 

Sncpiq) = S « Z! e(C.-)(C.-)p(9)> 

( 21 ) 

where a denotes the spin index. Similarly, the spin fluc¬ 
tuation Ssp{q) is defined by 



^ ^ ^ '^o’0’'6(C,o-)(Ccr')p(9)i (22) 

Qcrcr' 

where cr denotes the Pauli matrix. 

Finally, let us comment on the way to classify the fluc¬ 
tuation modes into charge and spin components. We clas¬ 
sify the eigenmodes according to the spin dependence as 
follows. The charge fluctuation Sn satisfies the relations 

( 5 (c|c.,-) = ( 5 (c|c^) 0 , ( 5 (c|cj^) = ( 5 (c|c.,-) = 0 . ( 23 ) 

On the other hand, the z component of the spin fluctua¬ 
tions satisfies 

( 5 (c|c.,-) =-,5 (c|c4,) 7^ 0 , , 5 (c|c^) = , 5 (c|c.,-) = 0 , ( 24 ) 

and the x and y components satisfy 

5 (c|ct) = i 5 (c|c 4 ,) = 0 , ( 5 (c|c^) = ±( 5 {c|c^)* 7^ 0 . ( 25 ) 

In addition, we categorize the eigenmodes into acoustic 
and optical ones according to the sublattice dependence. 
The acoustic mode satisfies the relation 

<5(41 c^i) = = (5(c4c^,)> (26) 

where rp denotes the position of the sublattice p. On the 
other hand, the optical mode satisfies the relation 

E 4 <c. 4 = 0 . ( 27 ) 

p=i 


III. RESULT 

In this section, we present the results for AIV2O4 and 
LiV204 obtained by the method in Sec. [Ill In Sec. IIII Al 
we show the LDA band structures and the tight-binding 
parameters. In Sec. IIII B[ we show the generalized sus¬ 
ceptibilities with the identification of the dominant fluc¬ 
tuations enhanced by electron correlations. 


A. LDA Band structure and MLWFs 

1 . Band structure 

We begin with the LDA calculations of the electronic 
band structures. We adopt a face-centered cubic primi¬ 
tive unit cell containing four V ions forming a pyrochlore 


lattice included in AIV2O4 and LiV2 04. The number of 
k points sampled is 4 x 4 x 4 . MLWFs are also composed 
by using the QMAS code. 

The left panels of Fig. [ 1 ] show the LDA band struc¬ 
tures for AIV2O4 and LiV2 04. The results are obtained 
for the experimental lattice structures at T = 800 K for 
AIV2O42S and 4 K for IjW 20 /^, respectively. The result 
for LiV2 04 well agrees with that in the previous study^S. 
In both cases, the relevant bands near the Fermi level 
are composed of the Sd t2g orbitals of V cations, ener¬ 
getically separated from other bands. We then perform 
the MLWF fitting for the t2g bands^i^. As shown in 
Fig. [ 1 ] the obtained MLWFs well reproduce the LDA 
band structure. From the MLWFs, we construct the non¬ 
interacting part of the multiband Hubbard Hamiltonian 
for each compound, which includes the electron transfer, 
trigonal crystal field splitting, and spin-orbit coupling. 


(a) 


> 




G 
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(b) 


> 
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FIG. 1 . (color online). Electronic states for (a) AIV2O4 and 
(b) LiV 2 04 . The Fermi level is set to zero. In the left pan¬ 
els, the curves denote the band structures obtained by the 
relativistic LDA calculations, while the results of the MLWF 
fitting are shown by the crosses. The right panels represent 
the DOS calculated by the MLWFs. The results decomposed 
into the aig and eg components are also shown. The inset in 
the right panel of (b) represents the first Brillouin zone for 
the pyrochlore lattice structure composed of V cations. 


The density of states (DOS) is shown in the right pan¬ 
els of Fig. [T] The results for AIV2O4 and LiV204 show 
contrasting behavior. In AIV2O4, the DOS is featureless 
near the Fermi level, and the oig and two Cg orbitals are 
almost equally occupied. On the other hand, in LiV204, 
the DOS near the Fermi level strongly depends on the 
energy, which is dominated by the Oig component. The 
difference of the DOS mainly comes from the difference 
in the transfer integrals. 
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2. Tight-binding parameter 

Tables U and mi show the transfer integrals calcu¬ 
lated by the MLWFs. For both compounds, the nearest- 
neighbor transfer is larger than that of the second and 
third neighbors. It is also common that the most dom¬ 
inant transfer is the dda component between nearest- 
neighbor sites, which is schematically shown in Fig. [ 2 ] 
However, the ratio of the diagonal dda component to 
the off-diagonal ones is largely different between the two 
compounds: the dda transfer is more than 30 times as 
large as the off-diagonal dxy-dy^ one in AIV2O4, whereas 
it is about twice of the dzx-dyz one in LiV204. 


For AIV2O4, as shown in Table IH^d), both of the off- 
diagonal components are nearly three times as large as 
the diagonal oig-aig one in magnitude, and comparable 
to the diagonal e^-Cg ones. In contrast, for LiV204, the 
diagonal Gg-Cg components are dominant, as shown in Ta- 
bleimd). The quantitative difference plays a decisive role 
in the contrasting behavior between the two compounds, 
as detailed below. 

Unlike the transfer integrals, the coupling constants 
of the one-body on-site terms are almost the same for 
the two compounds. For AIV2O4 and LiV204, we ob¬ 
tain the estimates the trigonal crystal field splitting 
Atrig = 53 meV and 61 meV, respectively, while spin- 
orbit coupling Asoi = 26 meV and 27 meV, respectively. 


TABLE I. Transfer integrals obtained from the 

MLWF analysis for AIV2O4 between (a) nearest-neighbor, 
(b) second-neighbor, and (c) third-neighbor sites in the basis 
of {dxy,dyz,dzx)- (d) represents the nearest-neighbor trans¬ 
fer integrals in the basis of (aig, Cg^, eg^ 2 )- The displace¬ 
ment vector S is (a/4, a/4,0), {a/2, a/4, —a/4), {a/2,a/2,Q), 
(a/4, a/4, 0) for (a), (b), (c), and (d), respectively, a is the 
lattice constant for the cubic unit cell of the pyrochlore lattice. 
The unit is meV. The values less than 1 meV are omitted. 
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TABLE 11. Transfer integrals for LiV 2 04. The notations are 
the same as in Table H) 
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B. Generalized Susceptibility and Eigenmode 
Analysis 

To study the fluctuations of AIV2O4 and LiV2 04, gen¬ 
eralized susceptibility is calculated. After the computa¬ 
tion of the static bare susceptibility, electron correlation 
is included at the level of RPA. By the eigenmode anal¬ 
ysis of the generalized susceptibility, charge, spin, and 
orbital fluctuations are obtained. 


1. Bare susceptibility 


Based on the non-interacting Hamiltonian constructed 
from the MLWFs, we first calculate the static bare sus¬ 
ceptibility Xalp-a'P'p'i^)- matrix size is 144 x 144 : 
four sublattices (only diagonal) with three orbital and 
two spin components [4 x (3 x 2 )^ = 144 ]. In the calcula¬ 
tions, we replace the integral of k over the first Brillouin 
zone by the summation for 32 ^ k grid points, while we 
take the summation of ojk analytically. 



FIG. 3. (color online). Eigenvalues of calculated by 

using the tight-binding parameters estimated by the MLWFs 
for (a) AIV2O4 and (b) LiV2 04 . The calculations are done 
at T = 10 K. 


FIG. 2. (color online). Schematic picture of edge-sharing 
VOe octahedra. The black and red spheres denote V and O, 
respectively. A pair of d^y orbitals is drawn to show the dda 
overlap. 

The difference of the transfer integrals between the two 
compounds is also seen in the (aig, e^) representation. 


Figure [ 3 ] shows all the eigenvalues of x^°^(g) for 
AIV2O4 and LiV204. The results look very differently: 
the eigenmodes are rather entangled in the entire spec¬ 
trum for AIV2O4, whereas the result for LiV204 shows 
sixteen eigenmodes well separated from the other modes. 
Analyzing the eigenvectors of the sixteen eigenmodes, we 
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find that they are dominantly in the aig orbital sector, 
which is anticipated from the DOS in Fig. [Ijb). 


2 ,. RPA calculation and eigenmode analysis of AIV2O4 

Next, we calculate the generalized susceptibility by in¬ 
cluding the effect of electron correlations in a perturba¬ 
tive way at the level of RPA. First, let us discuss the 
results for AIV2O4. We investigate the generalized sus¬ 
ceptibility obtained by RPA, while changing U, 

Jh, and V. Among the parameters, we find that V in¬ 
duces particularly interesting behavior, which might be 
related with the heptamer formation in AIV2O4, as dis¬ 
cussed below. 


(b) i-o 
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i 0.8 

I 0.7 
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FIG. 4. (color online), (a) V dependence of the maximum 
eigenvalues of for AIV2O4. (b) V dependence of the 

eigenvalues at go = (tt/IB, 7r/16,7r/16), which is the smallest 
wave number along the F-L line in the present calculations. 
In (a) and (b), we take U — 300 meV and Jh = 30 meV, 
while the other parameters are the same as in Fig. Ha), (c) 
a-bonding type charge fluctuations obtained from the eigen¬ 
mode analysis. The vertical axis represents the fluctuation of 
the electron density in arbitrary units. The histogram repre¬ 
sents the density fluctuations decomposed into the d^y, dyz, 
and dzx orbitals. The numbers 1-4 in the horizontal axis de¬ 
note the sublattices, (d) Schematic visualization of the fluc¬ 
tuations in the mode 1 in (a). 


Figure shows the largest eigenvalues of 
which are considerably enhanced by increasing V. As 
shown in Fig. Hb), the enhancement occurs in particu¬ 
lar eigenmodes, whereas all the other modes are almost 
insensitive to V. To clarify the nature of these enhanced 
fluctuations by V, we analyze the eigenvectors of the 
three quasi-degenerate modes^^. We find that the dom¬ 
inant fluctuations are in the charge sector. Figure He) 
shows the fluctuations of local electron densities. In all 
the three modes, the density fluctuation at one sublat¬ 
tice has the opposite sign to the other three, and the net 
density fluctuation vanishes in the four-site tetrahedron. 
Note that, while the densities at the sublattices 1 , 2 , and 


3 are suppressed in the modes 1 , 2 , and 3 , respectively, 
one can make the mode 4 in which the sublattice 4 is 
suppressed by a linear combination of the modes 1 - 3 . 

Interestingly, we find that the density fluctuations in 
these modes are strongly orbital dependent, as shown 
in Fig. He)- The orbital dependence indicates that the 
dominant fluctuations occur through the dda orbital on 
each bond. For instance, in the mode 1 , the charge den¬ 
sity in the dyz orbital is dominantly transferred between 
the sites 1 and 2 on the yz plane, which is regarded 
as the charge fluctuation through the dda orbital. The 
bond- and orbital-dependent fluctuations are schemati¬ 
cally shown in Fig. Hd)- Thus, the fluctuations of the 
three modes sensitive to V are of cr-bonding type. The 
importance of such cr-bonding states in the heptamer 
was suggested in the previous experimental and theo¬ 
retical studies^i^. Hence, our results for the dominant 
charge fluctuations in cr-bonding orbitals suggest that 
the inter-site Coulomb repulsion plays a role in the self¬ 
organization of heptamers in AIV2O4. 

We note that the values of V in the present RPA calcu¬ 
lations are rather large: in reality, the bare value of V will 
be considerably smaller than U. Nevertheless, our finding 
might be relevant to the heptamer formation in AIV2O4 
due to the following reason. The structural change as¬ 
sociated with the heptamer formation clearly indicates 
the importance of the Peierls-type electron-phonon cou¬ 
pling. Such inter-site phonons are known to give rise to 
an effective repulsive interaction for electrons: indeed, 
the integration of phonon degrees of freedom leaves the 
effective V term, together with other inter-site interac- 
tions^^. We regard that such effects are included in the 
value of V in the RPA analysis, although the realistic 
estimate is left for future study. 


3 . RPA calculation and eigenmode analysis of Li V2 O4 

Next, we discuss the results of x^^"^ LiV204. We 

here focus on the effect of U, while V leads to a differ¬ 
ent charge fluctuation from AIV2O4 as mentioned below. 
Figure shows U dependence of the maximum eigen¬ 
values of x^^^- The maximum eigenvalues are enhanced 
by U and become more dispersive. 

We plot [/ dependence of the largest sixteen eigenval¬ 
ues of x^^^ in Fig. [S{b)- Note that all of them are of oig 
character as discussed above. We find that nine eigen¬ 
modes are largely enhanced by U. We also elucidate that 
all the nine eigenmodes consist of spin fluctuations. Fig- 
ureHc) shows the three of them, which are the spin fluc¬ 
tuations in the x component. In all the three modes, the 
spin fluctuation at one sublattice has an opposite sign 
to the other three; the net Ss^ vanishes in the four-site 
tetrahedron. Similar situations are found also in the y 
and z components. Namely, the spin fluctuations satisfy 
the relation X]p=i = 0, where Ssp = (Ssp, Ss^, Ssp. 
Hence, we call them the optical-type spin fluctuations, 
whose schematic visualization is shown in Fig. H'i)' 
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FIG. 5. (color online), (a) U dependence of the maximum 
eigenvalues of for LiV 204 . (b) U dependence of the 

largest sixteen eigenvalues at qq. OPT (AGO) represents opti¬ 
cal (acoustic: uniform in the unit cell) type fluctuation. In (a) 
and (b), the calculations are done at Jn/U = 0.1 and V — 0 , 
while the other parameters are taken to be the same as in 
Fig-EKb). (c) OPT spin fluctuations obtained from the eigen¬ 
mode analysis. The vertical axis represents the spin fluctua¬ 
tions of the X component in arbitrary units. The other com¬ 
ponents are qualitatively similar, (d) Schematic visualization 
of the OPT spin fluctuations, in which the net spin fluctua¬ 
tion vanishes in the four-site tetrahedron, (e) T dependence 
of the maximum eigenvalues of X^^^(<7) at U — 300 meV, 
Jh = 30 meV, and P = 0. 


the zero-sum condition is often met in frustrated spin 
systems, our results suggest that the dominant fluctu¬ 
ations in LiV2 04 appear in the oig spins under strong 
geometrical frustration. 

On the other hand, we find that V enhances a charge 
fluctuation in the Oig orbital, where the net Oig den¬ 
sity fluctuation vanishes in the four-site tetrahedron (not 
shown here). The result is in sharp contrast to the case 
of AIV2O4, in which cr-type fluctuation is enhanced by 
P. 

Figure [SKe) shows T dependence of the maximum 
eigenvalues of x^^^ for LiV204. The dominant insta¬ 
bility is always in the optical-type spin fluctuations in 
the aig orbital. Although the mode is almost q inde¬ 
pendent at high T > 500 K, it develops the q depen¬ 
dence below ~ 300 K. At lower T below 100 K, the peak 
between the F-L line develops, being the dominant in¬ 
stability. We note that a similar peak appeared in the 
previous theoretical studyi^, while the model and the 
RPA treatment were different from ours. The growth of 
the incommensurate peak is potentially related with that 
observed in the inelastic neutron scattering for LiV2 04 
below 80 K— Nevertheless, further studies beyond the 
current weak-coupling approach are necessary to address 
the heavy fermion behavior, as the recent NMR experi¬ 
ment suggests local moments in the oig orbital even at 


high T^. 


IV. DISCUSSION 

It is worth noting that the electronic structure ob¬ 
tained by LDA and MLWF analyses plays an important 
role in the dominant fluctuations enhanced by electron 
correlations. In AIV2O4, three t2g orbitals almost equally 
contribute to the DOS near the Fermi level as shown in 
Fig. [IJa), and furthermore, the dda components in the 
transfer integrals are dominant as shown in TablelH These 
are directly reflected in the dominant cr-type charge fluc¬ 
tuation in Fig. 01 On the other hand, in LiV204, the 
diagonal components represented in the (oig, Cg) basis 
dominate the transfer integrals, as shown in Table IIIFd). 
This is closely related with the dominant oig spin fluc¬ 
tuation in Fig. [51 Thus, for both compounds, the or¬ 
bital channel for the dominant fluctuations are deduced 
from the transfer integrals composed from the LDA band 
structures, although the actual types of fluctuations are 
clarified only after performing RPA. 

Finally, we discuss the pressure effect on LiV204 in 
the present scheme. We performed the LDA calculation 
by reducing the lattice constant by 2% from the original 
value, which approximately corresponds to the situation 
at the pressure 6.3 GPapi^. By the energy optimization 
with respect to the u parameter, which controls the trig¬ 
onal distortion, we find that the value of u decreases from 
0.262 to 0 . 259 . Although it approaches that for AIV2O4 
(u = 0 . 251 )^, the oig component remains dominant in 
the DOS. Accordingly, our calculations for the general¬ 
ized susceptibility indicate that the dominant instability 
is still in the aig optical-type spin fluctuations and does 
not qualitatively change from ambient pressure. The 
results imply that the pressure-induced metal-insulator 
transition in LiV204 is caused by a different mechanism 
from the cluster formation in AIV2O4, and hence, the 
associated structural change is also likely to be different 
from AIV2O4. 


V. SUMMARY 

In summary, for a systematic understanding of the con¬ 
trasting behavior in two mixed valence spinels, AIV2O4 
and LiV2 04, we have analyzed charge-spin-orbital fluc¬ 
tuations in the multiband Hubbard models, whose one- 
body part is constructed by LDA calculations and MLWF 
analysis. From the eigenmode analysis of the general¬ 
ized susceptibility obtained by RPA, we found that, in 
AIV2O4, the CT-bonding type charge-orbital fluctuations 
are enhanced by the intersite repulsion, which provides 
a clue for understanding of the seven-site cluster forma¬ 
tion. In contrast, for LiV204, the optical-type spin fluc¬ 
tuation in the oig orbital is strongly enhanced by U at an 
incommensurate wave number at low temperatures. We 
also discussed the pressure effect on LiV204 and deduced 





















that the mechanism of pressure-induced metal-insulator 
transition might be different from the cluster formation 
in AIV2O4. 

As the present analysis is based on the LDA band 
structures and the electron correlations are included at 
the level of RPA, it is difficult to discuss the strong cor¬ 
relation effects, for instance, the large mass enhancement 
at low temperature in LiV204. It is necessary to adopt 
more sophisticated methods that can treat the renormal¬ 
ization of the band structures and occupations of each 
orbitals. This is left for future study. 
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